function list=intersected_triangles(triangle, points, level,edge_length)


index = find(abs(points(:,2)-level) < edge_length*1.5);
[x,y]=find(


